Quantum mechanical formalism for biological evolution 



O 

(N 

> 
O 

in 



o 

a 

I 

^— > 

» 

o3 



c 

o 
o 



> 

CO 
(N 

in 



X 



Gincstra Bianconi 1 and Christoph Rahmede 2 

1 Department of Physics, Northeastern University, Boston, Massachusetts 02115 USA 
2 Department of Physics and Astronomy, University of Sussex, Brighton, BN1 9QH,UK 

We study the evolution of sexual and asexual populations in general fitness landscapes. We 
find deep relations between the mathematics of biological evolution and the formalism of quantum 
mechanics. We give the general structure of the evolution of populations which is in general an 
off-equilibrium process that can be expressed by path integrals over phylogenies. These phylogenies 
are sums of linear lineages for asexual populations. For sexual populations instead, each lineage is 
a tree of branching ratio two and the path integral describing the evolving population is given by 
a sum over these trees. Finally, we show that the Bose-Einstein and the Fermi-Dirac distributions 
describe the stationary state of biological populations in simple cases. 
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I. INTRODUCTION 

The intriguing relation between evolutionary dynam- 
ics and statistical mechanics has attracted the in- 
terest of great evolutionary theorists like Fisher [l[ or 
Kimura Q that have related their results to the second 
principle of thermodynamics or to the theory of gases. 
Interestingly, the relation between evolutionary theory 
and quantum statistical mechanics is emerging from a 
series of independent works (5l4l3j that show a class of 
phase transitions occurring in evolution of haploid popu- 
lations and other evolving complex systems described by 
a Bose-Einstein condensation. In haploid populations, 
this transition is the quasi-species phase transition [l4r - 
Il9| in which a finite fraction of an asexual population 
ends up having the same genotype if the selective pressure 
is over a critical value and the mutation rate is smaller 
than a critical value. Moreover, in a recent paper [2(| 
it was shown that a condensation transition in the Bose- 
Einstein universality class occurs also in the evolution of 
diploid sexual populations in presence of epistatic inter- 
actions. When this condensation occurs, a finite fraction 
of pairs of genetic loci in epistatic interactions is fixed. 

The deep relation between evolutionary theory and the 
formalism of quantum mechanics extends also to the dy- 
namical description of biological evolution. The quantum 
spin-chain formalism has been shown to solve models of 
asexual evolution [2l[. Moreover, two papers have re- 
cently highlighted in s pec ific cases (asexual evolution in 
mean- field landscapes 12211 and asexual evolution in fluc- 
tuating environment |23f ) the role of path integrals in 
describing the temporal evolution of populations. 

In this paper we propose a general theoretical frame- 
work that reconciles all the cited results and extends 
both to asexual and sexual populations defined over com- 
plex fitness landscapes. The analogies that are present 
at mathematical level between biological evolution and 
quantum mechanics are highlighted. In particular we 
explore the different nature of the path integral calcu- 
lation that solves for the dynamical evolution of complex 
populations. The difference that resides between asex- 



ual and sexual evolution is that these phylogenies are 
sums of linear lineages for asexual populations, instead 
for sexual populations they are sums over trees of ances- 
tors, each individual having two parents, four grandpar- 
ents and so on. Finally we show that quantum statistics 
(both Fermi-Dirac and Bose-Einstein statistics) charac- 
terize the steady state distribution of sexual and asexual 
populations in simple cases. 

The paper is organized as follows: In section II we 
present the equations of evolution of asexual populations 
and in section III we present the equations and the struc- 
ture of their solution for sexual populations. In each of 
these sections we first study the off-equilibrium evolution 
with non-overlapping generations, then we study the evo- 
lution with overlapping generations and finally we char- 
acterize the steady state distribution of the evolutionary 
dynamics in simple cases. Finally in section IV we give 
the conclusions. 



II. EVOLUTION OF ASEXUAL POPULATIONS 

The genome of an asexual organism is formed by a 
single copy of each chromosome. If we indicate by 
i = 1, . . . , N a genetic locus, a given genotype is deter- 
mined by the allelic states {er} = (ax, 02, . . . ,Cj, . . . <jn) 
at each genetic locus i. The allelic state <%; at each ge- 
netic locus i can take 4 values corresponding to adenine, 
thymine, cytosine, guanine, i.e. <7j = 1,2,3,4. The pop- 
ulation evolves under the drive of selection that favors al- 
lelic configurations corresponding to higher reproduction 
rate, and mutations that increase the genetic variation 
in the population. We assume that the reproductive rate 
W({cr}), also called Wright fitness, of a genotype {a} is 
given by 



W({a}) = e-^« ff » 



(1) 



where U({a}) is the Fisher fitness and /3 is the selective 
pressure. If /3 = every genotype has the same reproduc- 
tive rate. If f3 3> 1 the difference in reproductive rate of 
genotypes having different U({o-}) is strongly enhanced. 
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A. Non-overlapping generations 

In the literature, usually the case of non-overlapping 
generations is studied 0, Hi] . The equations for asexual 
evolution are written in terms of the probability P ({cr}) 
that at generation t an individual of the population has 
genotype {cr}. The equation for the evolutionary dynam- 
ics in non-overlapping generations is given by [a, EH, Hi] 



P t+1 (M) 



where the normalization constant Z t is given by 



(2) 



^ = E E Q(W}\W})e- mW}) P'iW}) • (3) 



The matrix Q({c}|{cr'}) in Eq. (|2|) represents the prob- 
ability of mutations between genotype {cr'} to genotypes 
{c} in the next generation. Therefore, if fi is the muta- 
tion rate we have 



Q(wi{^'})=n[( i -^-^)+ 



(4) 



The solution of this evolutionary dynamics is given by 
a path integral over the phylogenies of the population 
Hi HI]. In fact if we iterate the Eq. © we find the 
solution 



?*({'}) = y E E -En» +1 iWr)a«)p o (Wo), 



(5) 



{cr} t _i {<r}t- 2 {<t}o t =° 



where the normalization constant is given by 



z t = n^ = E E •••EIl«(Wr + iiw.) e -^ ({CTK) J p (Wo) 

{cr} t {er}t-i {(t} r=0 



(6) 



r=0 



r 



The path integral in Eq. ([5]) can be directly calculated if 
we solve the eigenvalue problem 



M 



Wl{<r'} 



mW}) *n(W})] = A„^„({cr}) (7) 



where the operator M^n^/} ac ts on a generic function 
/({cr'}) according to the rule 

M W | K} /({a'}) = E Q(WIKI)/(K}) • (8) 

The eigenfunctions ir n ({cr}) of the eigenvalue problem 
(O are chosen to be normalized. Let us decompose the 
probability distribution P t {a}) on the basis of the eigen- 
functions ir n ({a}) according to the following expression, 
i.e. 



P 4 ({cr}) = Erfn(iK(W). 



(9) 



Using we obtain the dynamical equation for the coef- 
ficients d n (t), i.e. 



d n (t) = d„(0)A* 



(10) 



Therefore, asymptotically in time the population is con- 
verging to the distribution corresponding to the maximal 
eigenvalue Ao = max„A„. The partition function Z t is 
then given by 



(11) 



B. Overlapping generations 

We now write the equation of biological evolution in 
continuous time. Therefore we assume that at each time 
there can be a birth or a death process. We assume 
that the birth process depends on the reproductive rate 
IU({cr}) = cxp[— f3U({a})] and that the death process is 
a random drift. If we define the probability P({a}, t) that 
at time t and individual has genome {cr}, the dynamic 
equation of evolution of P({cr}, t) is given by 



dP(W},t) 
dt 



M 



-pu(W'}) 



P({a'},t) 



-P({a},t) 
(12) 
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where the operator M^ji^/i. is defined in ((SJ and 
Q({ct}|{(t'}) is given by (|3]). The partition function Z t in 
Eq. (|T2|) is given by 



z < = E E QCWIMJe-^"^^},*). (13) 

We observe here that the equation (fT2"j) can be reduced to 
the well known quasi-species equation (l5l - fl9| if we make 
a change of variables t —> if with dt/Z t = dt'. Let us now 
assume to know the solution of the eigenvalue problem 
(|7[l with A n describing the discrete spectrum of this prob- 
lem and TT n ({x}) the normalized eigenfunctions. If we 
decompose the function P({a},t) on the basis 7r„({<7}) 
of eigenfunctions , i.e. 



p(W,t) = 2c n (*Kn(M) , 



(14) 



the dynamical solution of Eq. (TT2"j) for the coefficients 
c n (t) is given by 



c(f) = exp [A„G(t) - t] c„(0) 



(15) 



In Eq. (|T5|) the function G(i) is defined through the 
function Z t according to the equation 



G(t) = / dt'— 
Jo &f 



(16) 



Using the definition for Z t given by Eq. (TIB")) together 
with Eqs. (|4]) and (|T5j) . we can close the self-consistent 
equations and uniquely determine the evolutionary dy- 
namics of the population. Therefore the partition func- 
tion Z t is given by 



(17) 



where the average (•) is performed over the functions 
c„(i) given by (fT5"]) , Finally, using (fT5"j) we can derive 
the equation obeyed by the partition function Z t = (A), 



1 d(\) 2 
2~dT~ 



«A 2 )-<A} 2 )>0. 



(18) 



This equation is related to the Fisher theorem of natural 
selection [l|, [23| and describes the fact that evolution is 
an off-equilibrium process. In fact, since (A) = {W'({cr})) 
is the average reproductive rate of the population, Eq. 
(|18p expresses the fact that this average reproductive 
rate always increases asymptotically in time as long as 
the environment does not change, i.e. the Fisher fitness 
function U{{a}) and the selective pressure remain con- 
stant in time. Eq. (|1 2|) shows surprising similarities with 
quantum mechanics that have been highlighted in a re- 
cent paper [25| presenting a unified framework between 
this evolutionary equation and stochastic quantization. 



C. The Bose-Einstein condensation in the 
Kingman model 

The Kingman model 0, Q is one of the most interest- 
ing stylized models of asexual evolution where the quasi- 
species transition is observed. In the framework of this 
model the quasi-species transition can be exactly mapped 
to the Bose-Einstein condensation in a Bosc gas. One in- 
teresting aspect of the Kingman model is that the evolu- 
tionary dynamics reaches an equilibrium due to the con- 
stant drive of random mutations. In the Kingman model 
each individual is assigned a single real parameter e > 
determining its reproductive rate W(e), i.e. 



W{e) 



-Pe 



(19) 



Moreover, in this model, when a mutation occurs a new 
offspring is generated with random fitness e drawn from a 
given distribution p(e). Therefore, instead of writing Eq. 
(|T2l for the probability P({a}) that an individual of the 
population has genotype {a}, we can write the equation 
for the probability V(e,t\to) that a random individual in 
the population is associated with a given Fisher fitness 
e at time t and had the last mutation at time to. The 
evolution of of P(e,t\to) is given by 



dP(e,t\t Q ) e-^'P(e,t\t ) 



P(e,t\to) . (20) 



dt Z t 
The partition function Zt in Eq. (|20|) is given by 

Z * = J dt ° J dee~ Pe P{e,t\t ) . 



(21) 



Therefore, the probability P(e,t\to) that an individual 
has fitness e at time t under the condition that the last 
mutation happened at time to, is given by the solution of 
©, i.e. 

P(e,t\to)^ap(e)e^o dtCKp[ - pe]{1 - a) ^' {t - to) . (22) 

Finally, integrating over to, we can evaluate the probabil- 
ity P(e, t) that an individual at time t has a given Fisher 
fitness e independently of to, i.e. 

P{e,t) - f dtoap{e)e^ dtc ^ [ -^ ]{1 - a) ^- (t - to \ 



In the case in which the distribution of the fitness has 
a finite support, one can easily find that he partition 
function Z t converges to a time independent value Z. 
Therefore the steady state solution for Pg(e) reached in 
the limit t — > oo is given by 



P* B (e)=pp(e) 



1 



1 



(23) 



where e~P^ B = (e~ ,3e ) /(1—fJ,). Therefore the probability 
that an individual has fitness e is determined by the Bose- 
Einstein distribution. The probability distribution Pg(e) 
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must satisfy the normalization condition 
1 = J deP* B {e)=p J 1 + J dep( 



(24) 



When p(e) is vanishing for e — > 0, the Bose-Einstein inte- 
gral in Eq. (|24|) can be limited from above. As a result, 
at high enough selection pressure and low enough muta- 
tion rate the system might undergo a condensation phase 
transition in the Bosc-Einstcin universality class. Below 
this phase transition a finite fraction of the individuals in 
the population shares the same genotype corresponding 
to the maximal fitness. This is one of the principal ex- 
amples that show the so called quasi-species transition: 
for low mutation rate and high selective pressure a fi- 
nite fraction of the population is found to have the same 
genotype. Interestingly a similar phase transition occurs 
also in evolving ecologies [l2| where the invasive species 
might strongly reduce the biodiversity, and in evolving 
models of complex networks 0, where there might be 
the emergence of super-hubs like Google in the World- 
Wide-Web. 



D. The Fermi-Dirac distribution in presence of 
negative selection 

The Fermi-Dirac distribution can be obtained in pres- 
ence of negative selection by a similar mechanism that 
generates the Bosc-Einstcin distribution in the King- 
man model. This model is mostly interesting because of 
the underlying symmetry between Fermi-Dirac and Bosc- 
Einstcin distribution. We assume that each individual is 
assigned a parameter e describing its adaptability to the 
environment. The birth probability is one, and new in- 
dividuals are generated at each time with parameter e 
drawn by a given distribution p(e). The death probabil- 
ity is given by a random drift with probability a and by 
a negative selection proportional to exp[/3e] with proba- 
bility 1 — a. Therefore the dynamics of evolution for the 
probability P(e,t\to) that in the population there is an 
individual with parameter e born at time to is given by 



dP(e,t\t Q ) 
dt 



(1 



~z7 



P(e,t\t ) 



The partition function Z± in Eq. (|25|) is given by 



Zt = Jo J deePtp ^ 1 ^ 



(25) 



(26) 



Solving Eq. (|25[) we find the probability that an individ- 
ual has fitness e at time t under the condition that it was 
born at time t , 

P(e,t\t ) = p ( e ) e /4 rft -PlMU-«)^-a( t - t0 ) _ (2?) 

The probability that an individual at time t has Fisher 
fitness e independently of to is 

P(e,t) = [\top(e)e f t° dtcM - M{1 - a] ^- a{t - to \ 



This system reaches steady state in the limit t — > oo with 
the equilibrium distribution given by the Fermi-Dirac dis- 
tribution Pp{e) 1 and in particular we have 



a eP^ fJ,F > + 1 



(28) 



with e^ F = (e^ e )a/(l - a). This distribution is the 
Fermi-Dirac distribution. The duality between Fermi- 
Dirac and Bose-Einstein distribution has been also recog- 
nized in the framework of evolving network models 0, Q 
and of models for evolving ecologies fl2j . 



III. EVOLUTION OF SEXUAL POPULATIONS 

In sexual populations, the somatic cells of diploid or- 
ganisms have two homologous copies of each chromo- 
some. In fact the genome of each diploid individual is 
given by the pairs of chromosomes A and B of the two 
haploid gametes coming from the father and from the 
mother of the individual. Let us suppose that each ga- 
mete is identified by N loci indicated with Latin letters 
i = 1, 2, . . . , N. If we indicate with cr^ the allelic state 
at each locus i, the gamete is characterized by the se- 
quence {a} = (<7i, 02, . . . , <7i, . . . , ctjv) with <ji = 1, 2, 3, 4 
indicating respectively the nucleotides adenine, thymine, 
cytosine, guanine. Given this description of the ga- 
metes, each individual is characterized by the sequence 
{^/} = [K,a?) . . . KVf ),...,«,<)} where 
af^ B indicates the allelic states in each parental ga- 
metes A/B. The reproductive rate, or Wright fitness, 
W({a A , <r B }) of a diploid individual depends on the pairs 
of chromosomes and can be expressed according to the 
relation 



W({a A 1 a B }) = e-r ju ^ A -° B » 



(29) 



where U({a A ,a B }) is the Fisher fitness of the individ- 
ual with genotype {a A ,a B } and (3 plays the role of the 
selective pressure. 

The gametic life-cycle describes the information trans- 
fer at next generations. Two gametes A/B generate a 
new individual {fertilization) and the new individual, if 
it reaches the reproductive state, carries the information 
and gives rise (by meiosis) to new gametes S±, S2, ■ ■ ■ S n 
with n = 0, 1, 2 . . . which continue the gametic life-cycle. 
The process of meiosis is a process of reduction in the ge- 
netic information of each diploid individual to generate 
gametes which have only half of the number of chromo- 
somes. During meiosis a process of recombination can 
occur with small probability at given locations (recombi- 
nation hotspots) on the chromosome. When a recombina- 
tion event occurs, homologous sites on two chromosomes 
can mesh with one another and may exchange genetic 
information. 

The evolution of diploid populations can be studied by 
similar techniques as used for the evolution of asexual 
haploid populations. In particular we might distinguish 
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between the the dynamics assuming non-overlapping or 
overlapping generations. 



ability P*({<j}) that a gamete has allelic state {a} at 
generation t. The evolution of this probability follows 
the equation 



A. Non-overlapping generations 

In the case in which the generations do not overlap we 
write the evolutionary dynamics in terms of the prob- 

I 



pt+ l ({*}) = M W}{W A i(rB} 



(30) 



This equation expresses the evolutionary principles in 
such a way that evolution favors on one hand the repro- 
duction of high fitness individuals and on the other hand 
recombination and mutation processes that enhance the 
variation in the population. These principles are the ori- 
gin of the stochastic nature of evolution. We note here 
that Eq. ([3"U|) is significantly different from Eq. ©. 
In fact Eq. (|30|) is a quadratic equation in P t ({cr}) 



whenever Eq. ([2]) is a linear equation. The operator 
Mr^iir^A^Bi introduced in Eq. (J3UJ) indicates the free 
recombination of genetic material and the mutations oc- 
curring when a new gamete is generated. In particu- 
lar the operator M.^ a y^ a A a B-j is defined as the average 
over the probability of free recombination and mutations 
Q({a}\{<T A ,<T B }). Therefore the action of the operator 
over a generic function f({a A ,a B }) is given by 



J 



M W | {ct a ct b } [/({ctV 3 })] 



E 



Q({a}\{a A ,a B })f({a A ,a B }) 



(31) 



r 



where 



Q({a}\{a A ,a B }) 

N r 

n 



1 — <-/ As 1 



-6(a u a B ) + ^n(a) 



(32) 



The functional expression of Q({a}\{a A ,a B }) given by 
Eq. (f3"2")l describe the fact that in sexual populations 
both, recombination processes and mutations (occurring 
with rate n) enhance the variation in the population. The 
partition function Z t in Eq. (|30p is given by 



Z t = 



£ e-^ u ^ A ^P({a A })P{{a B }) . (33) 

{<^},{a B } 

Similarly to what happens for asexual evolution, the so- 
lution of Eq. (|3T)|) is given by a path integral over the 
phylogenies of the diploid sexual populations. Neverthe- 
less, the phylogenetic histories associated to one genome 
are much more complex than in the asexual case. In 



asexual populations, every individual of the population 
is the descendent from a linear lineage. In fact every in- 
dividual has a single progenitor at any given generation. 
Instead the phylogenetic history of a single individual in 
a sexual population is by itself a tree of two parents, four 
grandparents and so on back in time. 

In Fig. [T]wc show the phylogenies of asexual organisms 
as linear lineages that combined together form a tree. In- 
stead in Fig. [5] we show a single phylogenetic history of a 
diploid individual which is by itself a tree. Note that for 
ancient generations the progenitors at the different leafs 
of the tree might correspond to the same individual. We 
indicate with Tit) the phylogenetic tree that starts at 
generation C?(0) and ends at generation Git)- Each gen- 
eration Q(i) is formed by 2 t_i individuals and the entire 
tree Tit) is formed by J\f(t) = (2* +1 — 1) individuals. In 
the following we indicate by Greek letters a, 7 the ga- 
metes, represented in the phylogenies represented in Fig. 
[2] by the links of the tree. Using this nomenclature we 
can write the path integral solving Eq. (|30l) . i.e. 



6 



pt ^ = l E E 



A/"(t) 



MeS(t-i) { CT } e e(t-2) MgS(o) cer(t) 



where the normalization constant is given by 



N(t) 



*= E E • 

{<r}eS(t) {<t}gG(*-i) {o-)ee(o) a£T(t) 



£ g({<r B }|{^^})e-W«^ ; P°(W 7 ). 



(35) 



7£S(0) 





FIG. 1: (Color online) The descendant of an individual in an 
asexual population can be described by a phylogenetic tree. 
The history that has given rise to each individual is neverthe- 
less a single path on this tree (indicated in the figure with red 
arrows). 



In Eqs. ([34])-([35]) we indicated with {c^crf } the geno- 
type of the individual that generates the gamete a. We 
can characterize the operator M.^ a y^ (J A (J By by studying 
how it acts on a base of normalized functions cf> n ({a}). 
In general we will have the relation 



FIG. 2: (Color online)The phylogenetic history for an indi- 
vidual in a sexually reproducing species is the tree of his an- 
cestors. It is interesting to observe that while the trees of 
ancestors realistically describes the evolutionary history that 
has given place to an individual, the global population of a 
species is not in general decreasing. In fact many of the pri- 
mordial ancestors of a given individual in a species are likely 
to coincide. 



Let us express the probability distribution P({a}) in the 
basis of the functions <fi n ({o'}) according to 



M 



{a}\{a A ,a B } 



(36) 



£(M) = £A.(t)tf»(M) 



(37) 
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Prom Eq. ([30]) we obtain that the system reaches a steady 
state when D n {t) = D* and Z t = Z with 



ZxD* = £a? s ZW 



(38) 



The detailed analysis of the implications of this solution 
in specific cases will be performed in subsequent publica- 
tions. 



B. Overlapping generations 

The probability P({cr}) that a gamete has an allelic 
configuration {cr}, in a continuous time evolution, satis- 
fies the dynamical equation 



gg(W) 

0i 



M 



MK<= 



e-P u ^ A '° B »P{{a A })P{{a B }) 



(W) 



where the average fitness is given by 



(39) 



Z t = ]T e- m{ ° A ' aB}) P(W A })P(W B }) • (40) 
{* A },{° B } 



The operator M^ji^a ct b} introduced in Eq. (|39[) in- 
dicates the free recombination of genetic material oc- 
curring when a new gamete is generated. In particular 
the operator M.^ (7 y^ rT A a B^ is defined as the average over 
the probability of free recombination Q({<r}\{a A ,cr B }). 
Therefore the action of the operator Mr,} \t a A tCr B\ over 
a generic function f({a A , cr B }) is given by Eq. (pi) with 
Q({^}\W A ^ B }) g ivcn by Eq. {32]). The equation for 
sexual evolution is much more complex than the one 
for asexual evolution since it is non-linear. In order to 
solve it we evaluate, as we did in the model with non- 
overlapping generations, the operator Muir^ a B^ on 
a basis of normalized functions </>„({cr}). The operator 
M{o-}|{o- A ,o- B }i wm a °t on any general basis of functions 



derive the condition that C* and Z have to satisfy, 

z x C« = E^ ^rsC*Cs 
rs 



Although there are many similarities between evolution 
of sexual and asexual populations, we observe also rele- 
vant differences. In particular in sexual evolution there is 
no equivalent of Eq. (|18p . therefore the average reproduc- 
tive rate is not necessarily increasing in the population. 



C. Bose-Einstein condensation in sexual 
populations driven by epistatic interactions 

Epistatic interactions [lij between genetic loci are the 
non-additive terms to the Fisher fitness U({a A }, {cr B }). 
Understanding the implications of evolution in presence 
of epistatic interactions is one of the most promising fields 
in biological evolution since it might capture essential fea- 
tures for the analysis of the genetic origin of diseases. In 
a recent paper [20( the evolution of sexual populations 
in presence of pairwise epistatic interactions between ge- 
netic loci, and absence of mutations, was considered. In 
the model studied in [2(| the loci i = 1, 2 . . . , N are linked 
in an epistatic network formed by L links. The epistatic 
interactions between pairs of loci have a role in determin- 
ing the fitness function that can be written according to 
the expression 



4>n{{&}) according to Eq. (|36j) . Therefore the solution of 
the sexual evolution is given by 



P({cr}) = ^C„(^„({cr}) 



(41) 



with the coefficients C„ (t) satisfying the non- linear equa- 
tion 



dC n (t) = ErsK S Cr(t)C s (t) 

dt Z t 



C n (t) 



(42) 



The system reaches a steady state only if all the coeffi- 
cients C„ and the partition function Z t converge to time- 
independent values (C*,Z, respectively). Using (|4"2"j) wc 



W({a A ,a B }) 



n 

<ij> 



-/3 Uij (af 



(44) 



where the product is extended to all genetic loci < i,j > 
linked in the epistatic network. If the network is locally 
tree-like, the general structure of the solution to the evo- 
lutionary equation (|39j) is given by 



H 

En 

h=l <i,j> 



(45) 



where < i,j > indicates all the pairs of linked nodes 
present in the epistatic network. 

In [2(| the steady states of the evolutionary dynamics 
of the type as in Eq. (|45|) with H = 1 in presence of 
epistasis and absence of mutations have been fully char- 
acterized. These steady states are multiple and the distri- 
bution of pairs of allelic states is given by a Bose-Einstein 
distribution. A Bose-Einstein condensation can occur in 
which a finite fraction of pairs of genetic loci fixes. 



IV. CONCLUSIONS 

In this paper we have outlined the similarities between 
biological evolution of sexual and asexual populations 
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and the formalism of quantum mechanics. Moreover, we 
have shown that the dynamics of biological evolution of 
non-overlapping generations is expressed in terms of path 
integrals. We have characterized the solution of the equa- 
tions for evolving sexual and asexual populations finding 
interesting similarities and significant differences. Finally 
we have observed the emergence of quantum statistics in 
special cases where the population is at stationarity. 

In the book What is life? [2(| Erwin Schrodinger pro- 
posed that a deep relation might exist between biological 
evolution and quantum mechanics. Since then this idea 



has fascinated biologists and physicists p7l - l30| . Never- 
theless we still lack a firm scientific basis for this sugges- 
tive proposal. We believe that this paper, relating the 
formalism of biological evolution to the one of quantum 
mechanics will open new perspectives to investigate fur- 
ther the relation between biological evolution and quan- 
tum mechanics. In the future we plan to study evolution 
in presence of specific fitness functions including cpistasis 
between genetic loci and in adaptive networks in which 
the fitness function changes with time in a non stationary 
fashion, describing the evolution of species. 
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